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ABSTRACT 



We present the results of three-dimensional global magnetohydrodynamic 
(MHD) simulations of the Parker-shearing instability in a differentially rotating 
torus initially threaded by toroidal magnetic fields. An equilibrium model of 
magnetized torus is adopted as an initial condition. When /3 = Pgas/Pm&g ~ 1 
at the initial state, magnetic flux buoyantly escapes from the disk and creates 
loop-like structures similar to those in the solar corona. Inside the torus, growth 
of non-axisymmetric magneto-rotational (or Balbus & Hawley) instability 
generates magnetic turbulence. Magnetic field lines are tangled in small scale 
but in large scale they show low azimuthal wave number spiral structure. After 
several rotation period, the system oscillates around a state with (3 ~ 5. We 
found that magnetic pressure dominated (/3 < 1) filaments are created in the 
torus. The volume filling factor of the region where f3 < 0.3 is 2-10%. Magnetic 
energy release in such low-/? regions may lead to violent flaring activities in 
accretion disks and in galactic gas disks. 
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1. INTRODUCTION 

Magnetic fields in differentially rotating disks play essential roles in the angular 
momentum transport which enable the accretion and various activities such as X-ray flares 
and jet formation. Motivated by the Skylab observations of the solar corona, Galeev, 
Rosner, & Vaiana (1979) proposed a model of magnetically structured corona in accretion 
disks consisting of X-ray emitting magnetic loops. The magnetic loops can be created due 
to the buoyant rise of magnetic flux tubes (or flux sheet) from the interior of the accretion 
disk. Matsumoto et al. (1988) carried out two-dimensional magnetohydrodynamic (MHD) 
simulations of the Parker instability (Parker 1966) in nonuniform gravitational fields which 
mimic those in accretion disks. They showed that when (3 = P gas /P ma , g ~ 1, a plane-parallel 
disk deforms itself into evacuated undulating magnetic loops and dense blobs accumulated 
in the valley of magnetic field lines. The effects of rotation and shear flow, however, were 
not included in their simulations. 

Shibata, Tajima, & Matsumoto (1990) carried out two-dimensional MHD simulations 
of the Parker instability including the effects of shear flow and suggested that magnetic 
accretion disks fall into two types; gas pressure dominated (high-/?) solar type disks, and 
magnetic pressure dominated (low-/?) cataclysmic disks. In high-/? {(3 > 1) disks, magnetic 
flux escapes from the disk more efficiently as f3 decreases. 

Several authors (Hawley, Gammie & Balbus 1995 ; Matsumoto & Tajima 1995 ; 
Brandenburg et al. 1995; Stone et al. 1996) have reported the results of three-dimensional 
local MHD simulations of magnetized accretion disks by adopting a shearing box 
approximation (Hawley et al. 1995). In differentially rotating disks, magnetorotational 
instability (Balbus & Hawley 1991) couples with the Parker instability (Foglizzo & Tagger 
1995). Since Parker instability grows for long wave length perturbations along magnetic 
field lines, non-local effects may affect the stability and nonlinear evolution. 

Results of global 3D MHD simulations including vertical gravity were reported by 
Matsumoto (1999) and Hawley (1999). By adopting an initially gas pressure dominated 
torus threaded by toroidal magnetic fields, Matsumoto (1999) showed that magnetic 
energy is amplified exponentially owing to the growth of the Balbus & Hawley instability 
and that the system approaches a quasi-steady state with (3 ~ 10. Matsumoto (1999) 
assumed (3 = constant in the torus at the initial state. When f3o ~ 1, the deviation from 
magneto-rotational equilibrium introduces large amplitude perturbations. 

In this paper, we present the results of 3D MHD simulations starting from an 
equilibrium MHD torus threaded by initially equipartition strength (f3 ~ 1) toroidal 
magnetic fields. 
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2. MODELS AND NUMERICAL METHODS 

The basic equations we use are ideal MHD equations in cylindrical coordinate system 
(r,<f>,z). We assume that the gas is inviscid and evolves adiabatically. Since we neglect 
radiative cooling, our numerical simulations postulates that the disk is advection-dominated 
(see Kato, Fukue, & Mineshige 1998 for a review). 

The initial condition is an equilibrium model of an axisymmetric MHD torus threaded 
by toroidal magnetic fields (Okada et al. 1989). We assume that the torus is embedded 
in hot, isothermal, non-rotating, spherical coronal gas. For gravity, we use the Newtonian 
potential. We neglect the self-gravity of the gas. At the initial state the torus is assumed 
to have a constant specific angular momentum L. 

We assume polytropic equation of state P = Kp 1 where K is a constant, and 7 is the 
specific heat ratio. According to Okada et al. (1989), we assume 

< = 5 = H(i "' 2y ~ 1 (i) 

where va is the Alfven speed and H is a constant. For normalization we take the radius 
r where the rotation speed L/r equals the Keplerian velocity vko = (GM/r Y^ 2 as unit 
radius and set po — vko — 1 at r = r®. Using these normalizations, we can integrate the 
equation of motion into the potential form; 

* = ~\ + |a + ^i v s + 2(7^1) ^ = *° = constant ' ( 2 ) 

where v 2 is the square of the sound speed, R = (r 2 + z 2 ) 1 / 2 and \I>o = ^(r ,0). By using 
equation (2), we obtain the density distribution. 

= [ max[M/ + l/^-L 2 /(2r 2 ),0] ] 1/(7 - 1) 
P {K^/^-mi + Po'r 2 ^} J U 
where f3 Q = 2K/H is the plasma (5 at (r,z) = (r ,0). The parameters describing the 
structure of the MHD torus are 7, (3q, L and K. In this paper we report the results of 
simulations for parameters @ — 1, 7 = 5/3, L — 1, and K = 0.05. The density of the halo 
at R = r is taken to be Phaio/Po — 10" 3 . The unit field strength is B = pl^ 2 v K0 . 

We solve the ideal MHD equations in a cylindrical coordinate system by using a 
modified Lax-Wendroff scheme with artificial viscosity. We simulated only the upper half 
space (z > 0) and assumed that at the equatorial plane, p, v r , v^, B r , B^, and P are 
symmetric and v z and B z are antisymmetric. The outer boundaries at r = r max and at 
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z = z max are free boundaries where waves can transmit. The singularity at R = is 
treated by softening the gravitational potential near the gravitating center (R < 0.2r ). 
The number of grid points is (N r , N^, N z ) = (200,64,240). To initiate non-axisymmetric 
evolution, small amplitude, random perturbations are imposed at t — for azimuthal 
velocity. 

3. NUMERICAL RESULTS 

Figure la shows the initial condition. Color scale denotes the density distribution and 
red curves depict magnetic field lines. Figure lb shows density distribution and velocity 
vectors at t — 6.2to where to is the orbit time to — ^^o/vko- As the matter which lost 
angular momentum accretes to the central region, the MHD torus becomes disk-like. 
Velocity vectors indicate that matter flows out from the disk. 

After the non-axisymmetric Balbus & Hawley instability grows, the inner region of the 
torus becomes turbulent. The turbulent motions tangle magnetic field lines in small scale 
and create numerous current sheets (or current filaments). Figure lc shows the magnetic 
field lines projected onto the equatorial plane and the density distribution at z — 0. Note 
that Figure lc shows only the inner region where —2.5 < x/r < 2.5 and —2.5 < y/r < 2.5. 
In large scales, magnetic field lines and density distribution show low azimuthal wave 
number spiral structure. Figure Id shows that magnetic loop structures similar to those 
in the solar corona are created. The yellow surfaces show strongly magnetized regions 
where \B\ — 0.1B and red curves show magnetic field lines at t — 7.5t . The magnetic 
loops buoyantly rise from the disk due to the Parker instability. Numerical results indicate 
that small loops which newly appeared above the photosphere (emerging magnetic loops) 
develop into expanding coronal loops. We can observe magnetic loops elongated in the 
azimuthal direction and loops twisted by the rotation of the disk. 

Figure 2 shows the isosurfaces of plasma-/? at t — 6.2t . Orange surfaces show strongly 
magnetized regions where (3 = 0.1. Yellow and green surfaces show the region where 
(3=1 and (3 = 10, respectively. Inside the disk, strongly magnetized, low-/? filaments are 
created. As we shall show below, low-/? region occupies a small fraction of the total volume. 
Intermittent magnetic structures (filamentary strongly magnetized regions) similar to those 
in the solar photosphere develop in the disk. 

Figure 3a shows the spatial average of the magnetic energy, log((i? 2 / Sir) / (P)) (dashed 
curve), log(-B 2 /(87rP )) (dash-dotted curve) and log((— BrB^/ Ak) / P ) (solid curve) averaged 
in 0.7 < r/r < 1.3 and < z/r < 0.3. This figure shows that magnetic energy in the 
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equatorial region decreases within a few rotation period due to the buoyant escape of 
magnetic flux, and that averaged plasma (5 oscillates quasi-periodically around (5 ~ 5. 
The ratio of averaged Maxwell stress to initial equatorial pressure, a t h is a t h — 0.07 when 
6t < t < 12t . Figure 3b shows the time evolution of the accretion rate 

/•2tt r0.3ro 

M acc = / / pv r rdzdip (4) 
jo jo 

at r = 0.3r , and of the outflow rate 

M out = / / pv z rdrdip (5) 

JO JO 

at z = 3.0r . Accretion rate increases and fluctuates around a mean value. Numerical 
results also indicate quasi-periodic ejection of the disk material. Figure 3c shows the 
Poynting flux (E x B/4ir) z which passes through the plane at z = 1.0, 2.0, 3.0r . After 
a few rotation period, magnetic flux is convected from the equatorial region to the disk 
surface. The volume filling factor of the region where f3 < 0.3 is shown in Figure 3d. 
Strongly magnetized region where (3 < 0.3 occupies about 8% of the total volume after 
a few rotation period. Later, the filling factor decreases to 2%. After about 10 rotation 
period, magnetic flux is regenerated in the disk and filling factor shows the second peak 
when the averaged is maximum. Following this second peak, Poynting flux at z = 2r 

increases, which suggest that magnetic flux escapes from the disk. 



4. DISCUSSION 

We showed by 3D global MHD simulations that when a differentially rotating torus 
is threaded by equipartition strength toroidal magnetic fields, magnetic loops emerging 
from the disk continue to rise and form coronal magnetic loops similar to those in the 
solar corona. Inside the disk, magnetic turbulence drives dynamo action which maintains 
magnetic fields and keeps the disk in a quasi-steady state with f3 ~ 5. We successfully 
simulated the formation process of the solar type disk by 3D direct MHD simulations. 

Numerical results indicate that in differentially rotating disks magnetic field lines 
globally show spiral structure with low azimuthal wave number but locally fluctuating 
components create numerous current sheets. We expect that magnetic reconnection taking 
place in current sheets generate hot plasmas which emit hard X-rays. When the disk is 
optically thin, such reconnection events may be observed as large amplitude sporadic X-ray 
time variations characteristic of low-states in black hole candidates (Kawaguchi et al. 1999). 

We found that inside the torus, filamentary shaped, locally strongly magnetized, 
low-/5 regions appear. Even when (5 ~ 5 in average, low-/5 regions where (5 < 0.3 occupy 



- 6- 



2-10% of the total volume. The low-/? filaments are embedded in high-/? plasma. Such an 
intermittent structure is common in magnetized astrophysical plasmas. Numerical results 
indicate that low-/? filaments are re-generated after they once disappear. 

Although we assumed point gravity suitable for accretion disks, numerical results can 
qualitatively be applied to galactic gas disks. Our numerical results suggest that although 
/? > 1 in average, low-/? filaments exist in galactic gas disks. Magnetic reconnection 
taking place in such low-/? regions may heat the interstellar gas and create hot, X-ray 
emitting plasmas. Tanuma et al. (1999) proposed a model that magnetic reconnection in 
strongly magnetized (~ 30/zG) regions in our Galaxy creates hot plasma which emit 7KeV 
component of Galactic Ridge X-ray Emission (GRXE). The low-/? filaments can also confine 
the hot plasma and prevent it from escaping from the Galaxy. 

We thank Drs. K. Shibata, T. Tajima, S. Mineshige, T. Kawaguchi and K. Makishima 
for discussion. Numerical computations were carried out by VPP300/16R at NAOJ. This 
work is supported in part by the Grant of Japan Science and Technology Corporation 
and the Grant-in-Aid of the Ministry of Education, Science, Sports, and Culture, 
Japan(07640348, 10147105). 
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Fig. 1. — The density distribution and magnetic field lines in a model starting from 
equipartition strength (/3q = 1) toroidal field, (a) The initial state. Color scale shows 
the density distribution and red curves depict magnetic field lines, (b) Density distribution 
at t — 6.2to and velocity vectors in xz-plane. (c) Density distribution at z = (color scale) 
and magnetic field lines projected onto the equatorial plane (—2.5 < x/ro < 2.5,-2.5 < 
y/ r o < 2.5). (d) Isosurfaces of magnetic field strength \B/B \ = 0.1 (yellow) at t — 7.5t . 
The blue surface shows the slice in xy-plane at z — 1.65r . The red curves show magnetic 
field lines. 

Fig. 2. — The isosurface of plasma (3 (—2.5 < r/ro < 2.5, —2.5 < z/tq < 2.5) at t — 6.2to- 
The orange surfaces show the strongly magnetized region where f3 = 0.1. The yellow and 
green surfaces show the isosurface (3 = 1 and f3 — 10, respectively. 

Fig. 3. — (a) Time development of the mean magnetic energy averaged in the region where 
0.7 < r/r < 1.3 and < z/r < 0.3. The dashed curve, dash-dotted curve and solid curve 
show log(( J B 2 /87r)/(P)), log(£ 2 /(87rP )>, and log(- J B rJ B^/(47rP )), respectively, (b) Time 
development of accretion rate at r = 0.3ro and outflow rate at z — 3.0ro- (c) Poynting flux 
which go through the plane z = 1.0, 2.0, 3.0ro- (d) Volume filling factor of the region where 
(3 < 0.3 in 0.7 < r/r < 1.3 and < z/r < 0.7. 
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